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Abstract 

The scattering lengths and effective ranges that describe low-energy nucleon-nucleon scattering are 
calculated in the limit of SU(3)-flavor symmetry at the physical strange-quark mass with Lattice 
Quantum Chromodynamics. The calculations are performed with an isotropic clover discretization 
of the quark action in three volumes with spatial extents of L ~ 3.4 fm, 4.5 fm and 6.7 fm, and with 
a lattice spacing of b ~ 0.145 fm. With determinations of the energies of the two-nucleon systems 
(both of which contain bound states at these up and down quark masses) at rest and moving in the 
lattice volume, Liischer's method is used to determine the low-energy phase shifts in each channel, 
from which the scattering length and effective range are obtained. The scattering parameters, in 
the % channel are found to be m,/*' = 9.50±g;^|±J;^ and m^r^ = 4.61±g;§f and in the 
3 5i channel are m^a^ = 7.45±gJ|±g;^ and m^r^ = 3.71±g |f ±g §|. These values are consistent 
with the two-nucleon system exhibiting Wigner's supermultiplet symmetry, which becomes exact 
in the limit of large- N c . In both spin channels, the phase shifts change sign at higher momentum, 
near the start of the t-channel cut, indicating that the nuclear interactions have a repulsive core 
even at the SU(3)-symmetric point. 
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I. INTRODUCTION 



Decades of experimental measurements of the nucleon-nucleon (NN) cross sections have 
resulted in precise phase-shifts and scattering parameters that provide the cornerstone of 
nuclear physics. These two-body interactions, when combined with multi-body interactions, 
dictate the low-energy spectrum and interactions of nuclei, and also the equation of state 
of nuclear matter at moderate densities. Further, these interactions are responsible for 
the fine-tunings that permeate nuclear physics, and are responsible for producing sufficient 
carbon in the universe to allow for the emergence of life. During the last several years, 
there has been a substantial effort to determine the NN interactions directly from quantum 
chromodynamics (QCD) using the numerical technique of Lattice QCD (LQCD) [1-14]. 
Steady progress is being made toward this objective, but calculations at the physical light- 
quark masses have not yet been performed, and essentially only one lattice spacing has 
been used in calculations. While calculations at the physical light-quark masses, that are 
extrapolated to the continuum limits and are performed in volumes (and temporal directions) 
that are much larger than the inverse Compton wavelength of the pion, are required to verify 
the LQCD technology and provide a rigorous underpinning of modern nuclear interactions, 
calculations at heavier pion masses are equally important in understanding and quantifying 
the fine-tunings in nuclear physics. It is crucial to understand how the fine-tunings in 
nuclear physics translate into constraints on the five relevant fundamental parameters in 
the standard model of particle physics: the three light-quark masses, and the strong and 
electromagnetic coupling constants, and the NN interaction provides the simplest place to 
begin this investigation. 

One reason that there are presently few LQCD calculations of NN interactions is the 
significantly greater complexity of multinucleon systems as compared with systems of sin- 
gle mesons and baryons. A second reason is that significant computational resources are 
required to generate high-quality ensembles of gauge field configurations at or near the 
physical light-quark masses, an effort that has only become practical with the availability 
of petascale computers, and, as yet, these ensembles are not at sufficiently large volume 
to be of use in nuclear physics. At heavier quark masses, the resources required to gen- 
erate ensembles of lattice gauge configurations and light-quark propagators are relatively 
small, the degradation of the signal-to-noise in multinucleon correlation functions is signifi- 
cantly reduced, and thermal effects are exponentially suppressed, compared to calculations 
at lighter pion masses. For these reasons, we performed LQCD calculations of a number 
of s-shell nuclei and hypernuclei with A < 5 at the SU(3)-flavor symmetry point with the 
physical strange quark mass giving m w = 805.9(0.6) (0.4) (8.9) MeV, at a lattice spacing of 
b ~ 0.145 fm, and in lattice volumes with spatial extents L ~ 3.4 fm, 4.5 fm and 6.7 fm [13]. 
In this work, we continue this study and explore the NN scattering phase shifts below the 
inelastic threshold and the associated scattering parameters relevant below the t-channel 
cut at the SU(3) symmetric point. 

II. OVERVIEW OF THE LATTICE QCD CALCULATIONS 

Three ensembles of isotropic gauge-field configurations, generated with a tadpole-improved 
Liischer-Weisz gauge action and a clover fermion action [15], are used in this work and have 
been used previously to calculate the lowest-lying levels of the s-shell nuclei and hypernu- 
clei [13]. This particular lattice-action setup follows closely the anisotropic clover action 
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TABLE I: Parameters of the ensembles of gauge-field configurations and of the measurements used 
in this work. The lattices have dimension L 3 x T, a lattice spacing b, and a bare quark mass b m q (in 
lattice units) generating a pion of mass m^. N SIC light-quark sources are used (as described in the 
text) to perform measurements on iV c fg configurations in each ensemble. The three uncertainties 
associated with the pion mass are statistical, fitting systematic and that associated with the lattice 
spacing, respectively. 



Label 


L/b T/b 


P 


b m q 


b [fm] 


L [fm] 


T [fm] 


77V [MeV] 


m-n-L 




iV cfg 


iV src 


A 


24 48 


6.1 


-0.2450 


0.145 


3.4 


6.7 


806.5(0.3)(0)(8.9) 


14.3 


28.5 


3822 


96 


B 


32 48 


6.1 


-0.2450 


0.145 


4.5 


6.7 


806.9(0.3)(0.5)(8.9) 


19.0 


28.5 


3050 


72 


C 


48 64 


6.1 


-0.2450 


0.145 


6.7 


9.0 


806.7(0.3)(0)(8.9) 


28.5 


38.0 


1905 


54 



of the ensembles generated by the JLab group that we have used in our previous calcula- 
tions [4, 11, 16-19]. The parameter tuning and scaling properties of this action will be dis- 
cussed elsewhere [20]. One level of stout smearing [21] with p = 0.125 and tadpole-improved 
tree-level clover coefficient csw = 1-2493 are used in the gauge- field generation. Stud- 
ies [20, 22, 23] of the partially-conserved axial-current (PCAC) relation in the Schrodinger 
functional indicate that this choice is consistent with vanishing 0(b) violations, leading to 
discretization effects that are essentially 0(b 2 ). The parameters of the ensembles are listed 
in Table I, and further details will be presented elsewhere [20]. As two-nucleon systems are 
the focus of this work, relatively large lattice volumes are employed for the calculations, with 
correspondingly large values of m w L and m n T. In order to convert the calculated energies 
from lattice units (l.u.) into physical units (MeV), a lattice spacing of b = 0.1453(16) fm has 
been determined for these ensembles of gauge-field configurations from the T spectrum [24] . 

The A^ c fg gauge configurations in each of the ensembles are separated by ten hybrid 
Monte-Carlo (HMC) evolution trajectories to reduce autocorrelations, and an average of 
A" src measurements are performed on each configuration. The quark propagators are con- 
structed with gauge-invariant Gaussian-smeared sources with stout-smeared gauge links. 
These sources are distributed over a grid, the center of which is randomly distributed within 
the lattice volume on each configuration, and the quark propagators are computed using the 
BiCGstab algorithm with a tolerance of 10~ 12 in double precision. Quark propagators, either 
unsmeared or smeared at the sink using the same parameters as used at the source, give rise 
to two sets of correlation functions for each combination of source and sink interpolating 
fields, labeled as SP and SS, respectively. The propagators are contracted to form nucleon 
blocks projected to fixed momentum at the sink for use in the calculation of the correlation 
functions. The blocks are denned as 

Bf (p,t;x ) =E eiP ' X ^ /l) ' i '( x ^^o)# ) ' / (x,i;xo)4 /3) ' fe '(x,i;x )4% , (1) 

x 

where is a quark propagator of flavor /, and the indices are combined spin-color indices 
running over % = 1, . . . , A^A",,. 1 The choice of the fa and the tensor depend on the 



To be specific, for a quark spin component i s = 1,...,N S and color component i c — 1,...,N C , the 
combined index % — N c (i s — 1) + i c . 
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spin and isospin of the nucleon under consideration. For our calculations we used the local 
interpolating fields constructed in Ref. [25], restricted to those that contain only upper spin 
components (in the Dirac spinor basis). This choice results in the simplest interpolating 
fields that also have the best overlap with the nucleon ground states. Blocks are constructed 
for all lattice momenta |p| 2 < 4 allowing for the study of multinucleon systems with zero or 
nonzero total momentum and with nontrivial spatial wave functions. Interpolating operators 
used at the sink have blocks with back-to-back momentum to access excited states. These 
interpolating operators have small overlaps onto the ground state, and it is found that the 
correlation functions from these back-to-back interpolating operators are as well fit by single 
states as are the ground-state correlation functions. More sophisticated methods such as 
Matrix Prony [26] and GPoF [27, 28] have been applied to the correlation functions and do 
not provide significant improvement in the extractions, indicating that the correlators with 
back-to-back momenta are close to orthogonal. Future investigations with a fully variational 
basis of operators are required to give us confidence in extracting states beyond the first 
excitation. The dispersion relations of the single mesons and nucleons on these ensembles 
have been examined in Ref. [13], and the infinite volume extrapolations of the masses are 
m n = 0.59426(12)(11) l.u. = 805.9(0.6)(0.4)(8.9) MeV and M N = 1.20359(41)(61) l.u. = 
1.635(0)(0)(18) GeV. 

III. NUCLEON-NUCLEON SCATTERING IN THE CHANNEL 

In contrast to the real world, there is a bound state in the ^o-channel at the SU(3)-symmetric 
point with a binding energy of B nn = 15.9(2.7) (2.7) (0.2) MeV [13] 2 . Two-nucleon correla- 
tion functions with total momentum |P| = 3 and |P| = 1 evaluated in three lattice volumes, 
with spatial extent L = 3.4 fm, 4.5 fm, 6.7 fm, are used to extract the binding energy of the 
di-neutron. While the results are found to be consistent, volume effects are observed in the 
smallest volume and the binding energy determined in the largest volume is taken as the 
infinite volume value. In the largest volume, the exponentially suppressed deviations from 
the infinite volume value are negligible for this binding energy. With one bound state in 
this channel, Levinson's theorem dictates that the phase shift is <5^ 5o ^(0) = n at threshold. 

While the location of the bound state does not correspond to a real value for the scattering 
phase shift, it does provide a real value of k cot 5, and a valuable constraint on the scattering 
parameters in the effective range expansion, which is valid below the t-channel cut starting at 
|k| = m n /2. Positively shifted energy eigenvalues of two-body systems in the lattice volumes 
can be related to scattering phase shifts using Luscher's relation [31-33] and its extension to 
boosted systems [34-37] provided they lie below the inelastic threshold. The effective mass 
plots (EMPs) associated with the two-nucleon energy, reduced by twice the nucleon mass, of 
the first excited states with |P| = are shown in fig. 1. Plateau's in the EMPs are found in all 
volumes, leading to clean extractions of the energy splitting, and also to values of the phase 



2 It should be noted that the HALQCD collaboration does not find a bound state in the channel nor 
3 5i — 3 Di coupled-channels at this pion mass [12]. Their results are arrived at through the solution of the 
Schrodinger equation with "potentials" calculated on the lattice, a method which is theoretically unsound 
in various ways [29, 30], and, moreover, is considerably less direct than determining binding energies using 
simple spectroscopy. 

3 Here, |P| denotes the magnitude of the total momentum in units of 2n/L. 
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FIG. 1: The EMPs of the first excited states with |P| = in the channel in the L = 3.4 fm, L = 
4.5 fm and L = 6.7 fm ensembles, respectively. Twice the nucleon mass has been subtracted from 
the energy. The dark (light) shaded regions correspond to the statistical uncertainty (statistical 
and systematic uncertainties combined in quadrature) of the fit to the plateau over the indicated 
time interval. 



shift via Luscher's eigenvalue equation at those energies in the L = 3.4 fm and L = 4.5 fm 
ensembles. The energy splitting in the L = 6.7 fm ensemble is consistent with zero, and 
therefore straddles the lowest singularity in the eigenvalue equation, thereby providing no 
meaningful constraint on the phase shift or k cot 5. The energy splittings are determined with 
a correlated ^-minimization fit of a constant to the plateau in the EMPs over the fit ranges 
shown in fig. 1. Jackknife resampling is used to generate the covariance matrix required 
to define x 2 from the correlation functions. Consistent results for the extracted values of 
the energies are obtained using both one- and two-state exponential fits to the correlation 
functions. In the figures, the inner shaded region corresponds to the statistical uncertainty 
in the fit, derived from the ^-minimization, and the outer uncertainty corresponds to the 
fitting systematic uncertainty combined in quadrature with the statistical uncertainty. The 
fitting systematic uncertainty is determined by varying the fitting interval over the extent of 
the plateau region and accommodating the range of central fit values. The EMPs associated 
with the first excited state with |P| = 1 are shown in fig. 2. The correlation functions in 




FIG. 2: The EMPs of the first excited state with |P| = 1 in the channel in the L = 3.4 fm, L = 
4.5 fm and L = 6.7 fm ensembles, respectively. Twice the nucleon mass has been subtracted from 
the energy. The dark (light) shaded regions correspond to the statistical uncertainty (statistical 
and systematic uncertainties combined in quadrature) of the fit to the plateau over the indicated 
time interval. 



the L = 4.5 fm and 6.7 fm ensembles are not sufficiently precise to provide a statistically 
meaningful phase shift, as the energies are near singularities in Luscher's eigenvalue equation. 
This demonstrates a significant difficulty in determining phase shifts in large volumes where 
the poles of the eigenvalue equation are collapsing to form the infinite-volume scattering 
continuum. Table II shows the results extracted from the lowest-lying continuum levels 
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TABLE II: Results from the lowest-lying continuum states in the So channel. 



Ensemble 


|P| 


bAE 


k \/m n 


k cot 5 j m n 


H°) 


24 3 x 48 
24 3 x 48 
32 3 x 48 



1 



0.0358(13) (16) 
0.1609(16) (37) 
0.0165(13) (22) 


0.3506(64)(78) 
0.7197(41)(93) 
0.2373(92) (96) 


n i7c; + .034 +0.043 
U - L 1 °-0.031-0.036 

n qn+0.07+0.15 
U - OU -0.07-0.17 

n mn+°- 031 + - 057 

U.UOU_ Q28-0 .046 


63.4(3.8)(4.7) 
-67(5)(11) 
83(7)(13) 



from each of the ensembles. 

Below the inelastic threshold, at |k| 2 = M^m^ + m^/4, where k is the magnitude of the 
three-momentum of each nucleon in the center-of-mass (CoM) frame, the s-wave scattering 
amplitude can be uniquely described by a single phase shift and more directly k cot 5. Near 
threshold, and more generally, below the t-channel cut, kcot5 has a power-series expansion 
in terms of the kinetic energy of the two-nucleons, 

kcot5 = — + *r|k| 2 + P|k| 4 + C(|k| 6 ) , (2) 

called the effective range expansion (ERE), where a is the scattering length (using the nuclear 
physics sign convention), r is the effective range and P is the shape parameter. While the 
range of possible values of the scattering length is unbounded, the size of the effective range 
and shape parameter are set by the range of the interaction. In fig. 3, the extracted values of 




(k/m,) 2 (kJm K f 

FIG. 3: k cot 5 in the ^So channel. The positive energy values are given in Table II and the negative 
energy value is determined from the di-neutron binding energy. The left panel is a two-parameter fit 
to the ERE, and the right panel is a three-parameter fit to the ERE, as described in the text. The 
inner (outer) shaded region corresponds to the statistical uncertainty (statistical and systematic 
uncertainties combined in quadrature). 

k cot5/m n given in Table II for |P| = and from the di-neutron binding energy are shown as 
a function of Ikp/m 2 . The three points shown in fig. 3 lie significantly below the t-channel 
cut and so the ERE of k cot 5 can be fit to define the phase shift throughout this kinematic 
regime. With three points to fit, two-parameter (left panel) and three-parameter (right 
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panel) fits to the ERE of kcot5/m n are performed and are shown as the shaded regions in 
fig. 3. 

The successful description by a two-parameter fit indicates small values of the terms that 
are higher order in the ERE, consistent with what is observed at the physical pion mass. 
The scattering length and effective range determined from the two-parameter fit are 

m n CSo) - q c()+0.78+1.10 „ r ( 1 S ) _ a fil +0.29+0.24 ( o\ 

rn^a — y.ou_ 69 _ 080 , rn^i — ^-OJ-_o.3i-o.26 > \°) 
corresponding to 

^So) _ O QQ+0.19+0.27 r (%) _ , i on+0.071+0.059 f (a\ 

a — ^•OO_ 017 _ 20 ™ , ' — i - i<3U -0.077-0.063 1111 ■ W 

The uncertainties associated with a^ 5 °) and r^ s< -^ are correlated, and their 68% confidence 
region is shown in fig. 4. The uncertainty in the scattering length is asymmetric as it is 

5.5, -n 



5.0 




4.0 



3 ' 5 7 8 9 10 VT 12 

m n a (,So) 

FIG. 4: The 68% confidence region associated with m^a^ 80 ^ and m v r^ s ^ in the channel. The 
inner region corresponds to statistical uncertainties and the outer region corresponds to statistical 
and systematic uncertainties combined in quadrature. 

the inverse scattering length that is the fit parameter. The shape parameter obtained from 
the three parameter fit to the ERE expansion is consistent with zero: Pm^. = — ll^Tg. The 
scattering length and effective range extracted from the three-parameter fit are consistent 
with the two-parameter fit, but with larger uncertainties. A full quantification of the the- 
oretical error in the determination of the ERE parameters requires more calculations than 
are currently available. 

The phase shift below the t-channel cut can be determined from these fit parameters, 
and is shown in fig. 5, along with the results of the LQCD calculations and the phase shift 
at the physical values of the quark masses. We expect the phase shift predicted by the ERE 
to deviate significantly from the true phase shift near the start of the t-channel cut, and this 
is indeed suggested by fig. 5. Like the phase shift at the physical point, the phase shift at 
the SU(3) symmetric point is found to change sign at larger momenta, consistent with the 
presence of a repulsive hard core in the NN interaction. 
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FIG. 5: The phase shift in the So channel. The left panel is a two-parameter fit to the ERE, while 
the right panel is a three-parameter fit to the ERE, as described in the text. The inner (outer) 
shaded region corresponds to the statistical uncertainty (statistical and systematic uncertainties 
combined in quadrature) in two- and three-parameter ERE fits to the results of the Lattice QCD 
calculation. The vertical (red) dashed line corresponds to the start of the t-channel cut and the 
upper limit of the range of validity of the ERE. The light (green) dashed line corresponds to the 
phase shift at the physical pion mass from the Nijmegen phase-shift analysis [38]. 



IV. NUCLEON-NUCLEON SCATTERING IN THE 3 Si - 3 D 1 COUPLED CHAN- 
NELS 

At the SU(3)-symmetric point, the deuteron is bound [13] by B = 19.5(3.6)(3.1)(0.2) MeV 
which, as with the bound di-neutron in the channel, provides a constraint on k cot 5. The 
deuteron has a d-wave component induced by the tensor interaction, however mixing between 
the s-wave and the d-wave is higher order in the ERE and first appears at the same order 
as the shape parameter [39]. Therefore, while the scattering length and effective range are 
purely s-wave, the shape parameter is contaminated by the d-wave admixture. The EMPs 
associated with the first excited states with |P| = are shown in fig. 6, and with |P| = 1 
are shown in fig. 7. The correlation functions calculated on the L = 6.7 fm ensemble have 
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FIG. 6: The EMPs of the first excited states with |P| = in the 3 Si channel in the L = 3.4 fm, L = 
4.5 fm and L = 6.7 fm ensembles, respectively. Twice the nucleon mass has been subtracted from 
the energy. The dark (light) shaded regions correspond to the statistical uncertainty (statistical 
and systematic uncertainties combined in quadrature) of the fit to the plateau over the indicated 
time interval. 

energies that are too close to, or straddle, the singularities of Liischer's eigenvalue equation 
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FIG. 7: The EMPs of the first excited states with |P| = 1 in the 3 Si channel in the L = 3.4 fm, L = 
4.5 fm and L = 6.7 fm ensembles, respectively. Twice the nucleon mass has been subtracted from 
the energy. The dark (light) shaded regions correspond to the statistical uncertainty (statistical 
and systematic uncertainties combined in quadrature) of the fit to the plateau over the indicated 
time interval. 



TABLE III: Results from the lowest-lying continuum states in the 3 Si channel. 



Ensemble 


|P| 


bAE 


\k\/m n 


24 3 x 48 





0.0306(16)(23) 


0.324(8)(12) 


24 3 x 48 


1 


0.142(23)(20) 


0.6708(60)(54) 


32 3 x 48 





0.0115(17)(23) 


0.198(15)(19) 


32 3 x 48 


1 


0.0788(24) (40) 


0.496(9)(14) 



k cot 5 1 



■^t H 78.6(4.7)(6.9) 



,+1.1+0.73 
'-4.0-9.3 



in C-+7.1+6.5 
" iz -°-6. 9-6.1 



-0.069(32) (43) 109(9)(13) 
-1(13)(22) 



and are not useful in determining the phase shift. The results extracted from the fits to the 
plateau regions in these EMPs are given in Table III. 
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FIG. 8: kcotS in the 3 Si channel. The left panel is a two-parameter fit to the ERE, and the right 
panel is a three-parameter fit to the ERE, as described in the text. The positive energy values 
are given in Table III and the negative energy value is determined from the deuteron binding 
energy. The inner (outer) shaded region corresponds to the statistical uncertainty (statistical and 
systematic uncertainties combined in quadrature). 



9 



In fig. 8, the extracted values of k cot 5/m n given in Table III and from the deuteron 
binding energy are shown as a function of |k| 2 /m 2 . Following the procedure used to analyze 
the results in the ^-channel, again with three points to fit, two-parameter (left panel) and 
three-parameter (right panel) fits to the ERE of kcotS/m^ are performed and shown as 
the shaded regions in fig. 8. The scattering length and effective range determined from the 
two-parameter fit are 

m n ( 3 Si) _ 7 4C -+0.57+0.71 ( 3 5i) _ o 7 i +0.28+0.28 (r\ 

III^U — I • z ±0_o.53-0.49 i ~ °- ' 1 -0.31-0.35 ) \°) 

corresponding to 

a™ = 1.82±g : lS±8:S fm , = 0.906^1°;^ fin , (6) 

and fig. 9 shows the 68% confidence region for the extracted values of a^ Sl ^ and r(' iSl \ The 
shape parameter obtained from the three parameter fit to the ERE expansion is consistent 
with zero: Pm 3 = 2lgig- Again the scattering length and effective range extracted from the 
three-parameter fit are consistent with the two-parameter fit, but with larger uncertainties. 
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FIG. 9: The 68% confidence region associated with m n a^ Sl ^ and m^r^ 31 ^ in the 3 S\ channel. The 
inner region corresponds to statistical uncertainties and the outer region corresponds to statistical 
and systematic uncertainties combined in quadrature. 



The phase shift below the t-channel cut can be determined from these fit parameters, and 
is shown in fig. 10, along with the results of the LQCD calculations and the phase shift at 
the physical point. As in the channel, the phase shift predicted by the ERE is expected 
to deviate significantly from the true phase shift near the t-channel cut, and this is seen in 
fig. 10. Like the 3 <Si phase shift at the physical point, and the phase shift we have obtained 
in the channel, the phase shift at the SU(3) symmetric point is found to change sign at 
larger momenta, consistent with the presence of a repulsive hard core in the NN interaction. 
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FIG. 10: The phase shift in the 3 S% channel. The left panel is a two-parameter fit to the ERE, while 
the right panel is a three-parameter fit to the ERE, as described in the text. The inner (outer) 
shaded region corresponds to the statistical uncertainty (statistical and systematic uncertainties 
combined in quadrature) in two- and three-parameter ERE fit to the results of the Lattice QCD 
calculation. The vertical (red) dashed line corresponds to the start of the t-channel cut and the 
upper limit of the range of validity of the ERE. The light (green) dashed line corresponds to the 
phase shift at the physical pion mass from the Nijmegen phase-shift analysis [38]. 



V. NUCLEON-NUCLEON EFFECTIVE RANGES 



Unlike the scattering length, the size of the effective range and the higher-order contributions 
to the ERE are set by the range of the interaction. The leading estimate of the effective range 
for light quarks is r ~ l/m 7T , and higher order contributions are expected to be suppressed 
by further powers of the light-quark masses. It is natural to consider an expansion of the 
product m^r in the light-quark masses. While the most general form of the expansion 
contains terms that are non-analytic in the pion mass [40-43], for instance of the form 
m q \ogm q , with determinations at only two pion masses (including the experimental value) 
a polynomial fit function is chosen, 

m n r = A + B + ... . (7) 

In fig. 11, the results of our LQCD calculations of m^r are shown, along with the experi- 
mental value in each channel and a fit to the form given in eq. (7). While the uncertainties 
in the lattice determinations are somewhat large compared to those of the experimental de- 
termination, it appears that there is modest dependence upon the light-quark masses. The 
fit values are 

4(^0) _ 1 q/io+0.080+0.079 r>CSo) _ A 00+0.55+0.59 p„ V -l 
— J--<J^O_o.080-0.083 > D — ^- zo -0. 56-0. 57 ^ e V 

A<*> = 0.726±g:8g±S:8S , = 3.70igf 7 j£g GeV- 1 . (8) 

The two-parameter fit is clearly over simplistic and more precise LQCD calculations are 
required at smaller light-quark masses to better constrain the light-quark mass dependence 
of the effective ranges. 
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FIG. 11: The NN effective range in the channel (left panel) and the 3 5i channel (right panel). 
The inner (outer) shaded region corresponds to the statistical uncertainty (statistical and systematic 
uncertainties combined in quadrature) in a two-parameter fit to the results of the Lattice QCD 
calculation and the experimental value. 



VI. FINE TUNINGS AND SU(4) SPIN-FLAVOR SYMMETRY 

At the physical values of the quark masses, the deuteron is an interesting system as it is 
much larger than the range of the nuclear force. Its binding energy is determined by the pole 
in the scattering amplitude in the 3 S\ — 3 D\ coupled channels. It is known very precisely at 
the physical light-quark masses, Bd = 2.224644(34) MeV, and recently LQCD calculations of 
the deuteron binding have been performed at unphysical light-quark masses [8, 11, 13, 14]. 
Given that both the scattering lengths and effective ranges calculated in this work are 
large compared with the pion Compton wavelength (which naively dictates the range of the 
interaction for light pions), we explore the naturalness of the two-nucleon systems. In this 
context, naturalness is defined by the length scales of the system as compared to the range 
of the interaction. By contrast, a fine-tuned quantity is one in which the length scales of 
the system are unnatural over a small range of parameters of the underlying theory. 
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FIG. 12: The left panel shows the ratio of the scattering length to effective range in the 3 <Si 
channel. The right panel shows the normalized deuteron binding momentum versus the pion 
mass [8, 11, 13, 14]. The black point denotes the experimental value. 
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The left panel of fig. 12 gives the ratio of the scattering length to effective range in the 
3 jSi channel as a function of the pion mass. As the effective range is a measure of the 
range of the interaction, this figure reveals that the deuteron is becoming more natural at 
heavier light-quark masses. In the right panel of fig. 12, the deuteron binding momentum 
7d (related to the binding energy by B d = ■yj/M N ) normalized to the pion mass is shown 
as a function of the pion mass. In the chiral regime one would expect that that jd scales 
as as suggested by effective field theory [44-51]. However, at the heavy up and down 
quark masses used here, naive expectations based on the uncertainty principle suggest that 
the deuteron binding momentum, if natural, would scale roughly as the inverse of the range 
of the interaction. As the ratio of 7^ to /77-7I- clS db function of is not constant, but rather 
is falling, we conclude that pion exchange is no longer the only significant contribution to 
the long-range component of the nuclear force, consistent with the meson spectrum found 
at these quark masses. 

While more precise calculations at these quark masses are desirable, and LQCD cal- 
culations at other light-quark masses and at other lattice spacings are required to make 
definitive statements, the present calculations suggest that the deuteron remains unnatural 
over a large range of light-quark masses. This would imply that the unnaturalness of the 
deuteron binding energy at the physical point is a generic feature of QCD with three light 
quarks and does not result from a fine-tuning of their masses. If subsequently confirmed, 
this would be a very interesting result. 
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FIG. 13: The left panel shows the ratio of the scattering length to effective range in the 
channel. The right panel shows the normalized di-neutron binding momentum versus the pion 
mass [8, 11, 13, 14]. 



The ^So channel is unnatural at the physical point with a very large scattering length, 
but the system appears to be more natural at heavier pion masses. Nonetheless, as shown in 
fig. 13 (left panel), the scattering length is approximately twice the effective range at a pion 
mass of m n ~ 800 MeV, similar to the 3 <Si channel. In the right panel of fig. 13, the di-neutron 
binding momentum 7 nn (related to the binding energy by B nn = ■y^ in /M N ) normalized to the 
pion mass is shown as a function of the pion mass. As in the 3 Si channel, it appears that 
the pion is not providing the only significant contribution to the long-range component of 
the nuclear force. However, in contrast to the channel, the ^o-channel is clearly finely- 
tuned at the physical light-quark masses. The range of light-quark masses over which it is 
fine-tuned requires further LQCD calculations to determine, and eventual consideration of 
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isospin violating effects due to light-quark mass differences and electromagnetism. However, 
given the experimental determinations of the nn, np and pp scattering lengths, these effects 
are expected to be small. 

It is interesting to note that the ratio of the scattering length to the effective range in 
the two channels have very similar values at the quark masses used in this work: 

a^/r^ = 2.06«;^+0;25 , a (*o)/ r (fco) = 2 03+0.23+0.29 ? (Q) 

and that the scattering lengths in the two channels, and also the effective ranges, are within 
~ 20% of each other. In the large-N c limit of QCD, the nuclear forces in the two spin 
channels are equal up to corrections suppressed by 0(1/ N^) [52], and the two channels 
transform in the 6 of the Wigner SU(4) symmetry. In addition, inequalities for the binding 
energies of light nuclei in the Wigner-symmetry limit have been found in Ref. [53]. The 
closeness of the values of the scattering parameters at m 7r ~ 800 MeV is consistent with the 
expectations of the large-N c limit of QCD. 

VII. CONCLUSIONS AND DISCUSSIONS 

We have presented the results of Lattice QCD calculations of low-energy NN scattering 
phase-shifts and scattering parameters at the SU(3) symmetric point with a pion mass 
of ~ 800 MeV. For the first time, the effective ranges of the NN interactions have 
been determined using lattice QCD. The calculated scattering lengths and effective ranges 
indicate that the pion is not the dominant contribution to the long range part of the nuclear 
force at these large light-quark masses, as anticipated from the single-hadron spectrum. In 
both spin channels, the NN phase shifts change sign at higher momentum, near the start 
of the t-channel cut, indicating that the nuclear interactions have a repulsive core even 
for heavier quark masses. This suggests that the form of the nuclear interactions, and the 
effective potentials that will reproduce the scattering amplitude below the inelastic threshold, 
is qualitatively similar to the phenomenological potentials that describe the experimental 
scattering data at the physical pion mass. 

Both spin channels are, in a sense, more natural at m n ~ 800 MeV, where both satisfy 
a/r ~ +2.0, than at the physical pion mass where a^ s °^ /A 1 ^ ~ —8.7 and a^ Sl ^ /r^ Sl ^ ~ 
+3.1. The relatively large size of the deuteron compared with the range of the nuclear forces 
may persist over a large range of light-quark masses, and therefore might, in fact, not be 
usefully regarded as a fine-tuning in n/ = 2 + 1 QCD, but rather a generic feature. The 
channel, in contrast, is finely tuned at the physical light-quark masses and it remains to be 
seen over what range of masses this persists. 

Our calculations were performed at a single pion mass with one lattice spacing and in the 
absence of electromagnetic interactions. It should be stressed that in the presence of fine- 
tuning, as in the channel at the physical point, lattice-spacing artifacts can be enhanced 
with respect to expectations based on naive dimensional analysis and scaling arguments. In 
order to fully explore the behavior of the scattering phase shifts and scattering parameters 
with fully quantified uncertainties, along with the issues of spin-flavor symmetry and fine 
tunings, calculations at multiple lattice spacings and smaller light-quark masses are essential 
and are planned for the future. 



14 



Acknowledgments 



We thank R. Edwards and B. Joo for help with QDP++ and Chroma [54]. We ac- 
knowledge computational support from the USQCD SciDAC project, the National Energy 
Research Scientific Computing Center (NERSC, Office of Science of the US DOE, DE-AC02- 
05CH11231), the UW HYAK facility, LLNL, the PRACE Research Infrastructure resource 
CURIE based in France at the Tres Grand Centre de Calcul, TGCC, and the NSF through 
XSEDE resources under grant number TG-MCA06N025. SRB and PJ were partially sup- 
ported by NSF continuing grant PHY1206498. In addition, SRB gratefully acknowledges 
the hospitality of HISKP and the Mercator programme of the Deutsche Forschungsgemein- 
schaft. The work of AP is supported by the contract FIS2011-24154 from MEC (Spain) and 
FEDER. H-WL and MJS were supported in part by the DOE grant DE-FG03-97ER4014, 
and the NSF MRI grant PHY-0922770 (HYAK). KO was supported in part by DOE grants 
DE-AC05-06OR23177 (JSA) and DE-FG02-04ER41302. WD was supported by the U.S. 
Department of Energy under cooperative research agreement Contract Number DE-FG02- 
94ER40818 and by DOE Outstanding Junior Investigator Award DE-SC000-1784. The work 
of TL was performed under the auspices of the U.S. Department of Energy by LLNL under 
Contract DE-AC52-07NA27344. The work of AWL was supported in part by the Director, 
Office of Energy Research, Office of High Energy and Nuclear Physics, Divisions of Nu- 
clear Physics, of the U.S. DOE under Contract No. DE-AC02-05CH11231. MJS thanks the 
Alexandar von Humboldt foundation for the award that enabled his visit to the University 
of Bonn, and the kind hospitality of Ulf MeiBner and the Helmhotz-Institut fur Strahlen- 
und Kernphysik at the University of Bonn. 



[1] M. Fukugita, Y. Kuramashi, H. Mino, M. Okawa and A. Ukawa, Phys. Rev. Lett. 73, 2176 

(1994) [arXiv:hep-lat/9407012]. 
[2] M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino and A. Ukawa, Phys. Rev. D 52, 3003 (1995) 

[arXiv:hep-lat/9501024] . 

[3] S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage, Phys. Rev. Lett. 97, 012001 (2006) 
[arXiv:hep-lat/0602010]. 

[4] S. R. Beane et al. [NPLQCD Collaboration], Phys. Rev. D 81, 054505 (2010) [arXiv:0912.4243 
[hep-lat]]. 

[5] N. Ishii, S. Aoki and T. Hatsuda, Phys. Rev. Lett. 99, 022001 (2007) [arXiv:nucl-th/0611096]. 
[6] S. Aoki, T. Hatsuda and N. Ishii, Comput. Sci. Dis. 1, 015009 (2008) [arXiv:0805.2462 [hep- 
ph]]. 

[7] S. Aoki, T. Hatsuda and N. Ishii, Prog. Theor. Phys. 123, 89 (2010) [arXiv:0909.5585 [hep- 
lat]]. 

[8] T. Yamazaki, Y. Kuramashi, A. Ukawa, arXiv: 1105. 1418 [hep-lat]. 

[9] T. Yamazaki, Y. Kuramashi, A. Ukawa, Phys. Rev. D 81, 111504 (2010) [arXiv:0912.1383 
[hep-lat]]. 

[10] P. de Forcrand and M. Fromm, Phys. Rev. Lett. 104, 112005 (2010) [arXiv:0907.1915 [hep- 
lat]]. 

[11] S. R. Beane et al. [NPLQCD Collaboration], Phys. Rev. D 85, 054511 (2012) [arXiv: 1109.2889 
[hep-lat]]. 



15 



T. Inoue et al. [HAL QCD Collaboration], Nucl. Phys. A 881, 28 (2012) [arXiv:1112.5926 
[hep-lat]]. 

S. R. Beane, E. Chang, S. D. Cohen, W. Detmold, H. W. Lin, T. C. Luu, K. Orginos and 

A. Parreno et al, arXiv:1206.5219 [hep-lat]. 

T. Yamazaki, K. -i. Ishikawa, Y. Kuramashi and A. Ukawa, Phys. Rev. D 86, 074514 (2012) 
[arXiv: 1207.4277 [hep-lat]]. 

B. Shcikholeslami and R. Wohlert, Nucl. Phys. B 259, 572 (1985). 

S. R. Beane et al, [NPLQCD Collaboration], Phys. Rev. D 80, 074501 (2009) [arXiv:0905.0466 
[hep-lat]]. 

S. R. Beane et al. [NPLQCD Collaboration], Phys. Rev. Lett. 106, 162001 (2011) 
[arXiv:1012.3812 [hep-lat]]. 

S. R. Beane et al, [NPLQCD Collaboration], arXiv: 1103.2821 [hep-lat]. 

S. R. Beane, E. Chang, S. D. Cohen, W. Detmold, H. -W. Lin, T. C. Luu, K. Orginos 
and A. Parreno et al, [NPLQCD Collaboration], Phys. Rev. Lett. 109, 172001 (2012) 
[arXiv: 1204.3606 [hep-lat]]. 

W. Detmold, R. Edwards, B. Joo, T. Luu, S. Meinel, K. Orginos, D. Richards and A. Walker- 
Loud, in preparation. 

C. Morningstar and M. J. Peardon, Phys. Rev. D 69, 054501 (2004) [hep-lat/0311018]. 

R. Hoffmann, A. Hasenfratz and S. Schaefer, PoS LAT 2007 (2007) 104 [arXiv:0710.0471 
[hep-lat]]. 

R. G. Edwards, B. Joo and H. -W. Lin, Phys. Rev. D 78, 054501 (2008) [arXiv: 0803. 3960 
[hep-lat]]. 

S. Meinel, private communication (2012). 

S. Basak et al. [Lattice Hadron Physics (LHPC) Collaboration], Phys. Rev. D 72, 074501 
(2005) [hep-lat/0508018]. 

S. R. Beane, W. Detmold, T. C. Luu, K. Orginos, A. Parreno, M. J. Savage, A. Torok and 
A. Walker-Loud, Phys. Rev. D 79, 114502 (2009) [arXiv:0903.2990 [hep-lat]]. 
C. Aubin, K. Orginos, [arXiv: 1010.0202 [hep-lat]]. 
K. Orginos PoS(Lattice 2010)118 

S. R. Beane, W. Detmold, K. Orginos and M. J. Savage, Prog. Part. Nucl. Phys. 66, 1 (2011) 

[arXiv: 1004.2935 [hep-lat]]. 

M. C. Birse, arXiv: 1208.4807 [nucl-th]. 

M. Liischer, Commun. Math. Phys. 105 (1986) 153. 

M. Liischer, Nucl. Phys. B 354, 531 (1991). 

S. R. Beane, P. F. Bedaque, A. Parreno and M. J. Savage, Phys. Lett. B 585, 106 (2004) 
[hep-lat/0312004]. 

K. Rummukainen and S. A. Gottlieb, Nucl. Phys. B 450, 397 (1995) [hep-lat/9503028]. 

C. h. Kim, C. T. Sachrajda and S. R. Sharpe, Nucl. Phys. B 727, 218 (2005) [hep-lat/0507006] . 
N. H. Christ, C. Kim and T. Yamazaki, Phys. Rev. D 72, 114506 (2005) [hep-lat/0507009]. 
Z. Davoudi and M. J. Savage, Phys. Rev. D 84, 114502 (2011) [arXiv:1108.5371 [hep-lat]]. 
Nijmegen Phase Shift Analysis, http://nn-online.org/. 

J. -W. Chen, G. Rupak and M. J. Savage, Nucl. Phys. A 653, 386 (1999) [nucl-th/9902056]. 
S. Weinberg, Phys. Lett. B 251, 288 (1990). 
S. Weinberg, Nucl. Phys. B 363, 3 (1991). 

D. B. Kaplan, M. J. Savage and M. B. Wise, Phys. Lett. B 424, 390 (1998) [nucl-th/9801034]. 
D. B. Kaplan, M. J. Savage and M. B. Wise, Nucl. Phys. B 534, 329 (1998) [nucl-th/9802075]. 



16 



[44] S. R. Beane, P. F. Bedaque, M. J. Savage, and U. van Kolck, Nucl. Phys. A 700, 377 (2002) 
[nucl-th/0104030]. 

[45] S. R. Beane and M. J. Savage, Nucl. Phys. A 713, 148 (2003). [hep-ph/0206113]. 
[46] E. Epelbaum, U.-G. MeiBner, W. Glockle, Nucl. Phys. A 714, 535 (2003) [nucl-th/0207089]. 
[47] S. R. Beane and M. J. Savage, Nucl. Phys. A 717, 91 (2003) [nucl-th/0208021]. 
[48] E. Braaten, H.-W. Hammer, Phys. Rev. Lett. 91, 102002 (2003) [nucl-th/0303038]. 
[49] V. V. Flambaum and R. B. Wiringa, Phys. Rev. C 76, 054002 (2007) [arXiv:0709.0077 [nucl- 
th]]. 

[50] J.-W. Chen et al, Phys. Rev. C 86, 054001 (2012) [arXiv: 1012.0453 [nucl-th]]. 

[51] M. E. Carrillo-Serrano, I. C. Cloet, K. Tsushima, A. W. Thomas and I. R. Afnan, Phys. Rev. 

C 87, 015801 (2013) [arXiv: 1208.3009 [nucl-th]]. 
[52] D. B. Kaplan and M. J. Savage, Phys. Lett. B 365, 244 (1996) [hep-ph/9509371]. 
[53] J. -W. Chen, D. Lee and T. Schafer, Phys. Rev. Lett. 93, 242302 (2004) [nucl-th/0408043]. 
[54] R. G. Edwards and B. Joo, Nucl. Phys. Proc. Suppl. 140 (2005) 832. 



17 



